Harmonic vs. subharmonic patterns in a spatially forced oscillating chemical reaction 
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The effects of a spatially periodic forcing on an oscillating chemical reaction as described by the 
Lengyel-Epstein model are investigated. We find a surprising competition between two oscillating 
patterns, where one is harmonic and the other subharmonic with respect to the spatially periodic 
forcing. The occurrence of a subharmonic pattern is remarkable as well as its preference up to rather 
large values of the modulation amplitude. For small modulation amplitudes we derive from the model 
system a generic equation for the envelope of the oscillating reaction that includes an additional 
forcing contribution, compared to the amplitude equations known from previous studies in other 
systems. The analysis of this amplitude equation allows the derivation of analytical expressions even 
for the forcing corrections to the threshold and to the oscillation frequency, which are in a wide range 
of parameters in good agreement with the numerical analysis of the complete reaction equations. 
In the nonlinear regime beyond threshold, the subharmonic solutions exist in a finite range of the 
control parameter that has been determined by solving the reaction equations numerically for various 
sets of parameters. 

PACS numbers: 82.40.Ck, 47.20.Ky, 47.54.-r 
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I. INTRODUCTION 

Studying the response of pattern-forming systems 
with respect to an external stimulus provides a pow- 
erful method to investigate the inherently nonlinear 
mechanism of self-organization in various systems under 
nonequilibrium conditions. Thermal convection fil and 
electroconvection in nematic liquid crystals 0, EST! are 
two systems, where the effects of spatially periodic forc- 
ing have been investigated rather early. Further on, the 
effects of forcing on stationary bifurcations have been 
studied extensively in many different systems [E IE IS 
El EH EJL EH EH - an d this branch of nonlinear science has 
also evolved to forcing studies on oscillatory media and 
traveling waves [IllllElElEllllllSllilllllllll. 

Recently, the response behavior of patterns with re- 
spect to a combination of spatial and temporal forcing 
has attracted a great deal of attention because of the de- 
velopment of flexible forcing techniques using iUumma- 
tion, as for instance in photosensitive chemical reactions 
[2^. I24L I25L 12a . |27L [2^1 or in electroconvection in nematic 
liquid crystals |23, y(| • For the photosensitive chlorine 
dioxide-iodine-malonic acid (CDIMA) reaction, as de- 
scribed by the so-called Lengyel-Epstein model [3ll l32| , 
one finds in a large parameter range Turing patterns. In 
particular, their response with respect to a forcing of a 
traveling wave type, which is spatially resonant or near- 
resonant with respect to the characteristic wavelength of 
the Turing pattern, exhibits a number of new phenom- 
ena and has therefore attracted considerable attention 
recently [Mill- 

The Lengyel-Epstein model also exhibits a spatially ho- 
mogeneous and supercritical Hopf bifurcation |3lLl33.l34| 
similar to the one found in other chemical reactions. In 
the present work, we investigate the response of the Hopf 
bifurcation of this model with respect to a spatially pe- 



riodic but time-independent illuminating forcing, which 
enters the Lengyel-Epstein model additively. Beyond the 
instability of the homogeneous chemical reaction, we find 
a surprising competition between a temporally oscillating 
and spatially modulated reaction that is harmonic with 
respect to the spatially periodic forcing and another one 
which is subharmonic. Besides its occurrence close to the 
threshold of the Hopf bifurcation, also the preference of 
the subharmonic pattern up to large forcing amplitudes 
is remarkable as well. 

From a technical point of view, our analysis is related 
to a previous study of a complex Ginzburg-Landau equa- 
tion corresponding to a spatially periodic modulation of 
a temporally resonant forcing of a chemical reaction |25| . 
While the forcing in the previous studies entered the 
model equation multiplicatively, the forcing enters the 
Lengyel-Epstein model additively. Close to the thresh- 
old of the Hopf bifurcation, we also reduced the Lengyel- 
Epstein model to a universal equation for the amplitude 
of the oscillations by using a multiple scale perturbation 
technique [§3 ■ We find that the forcing contribution 
occurs also multiplicatively in the resulting amplitude 
equation, but the forcing contribution to the complex 
amplitude equation has a different form, compared to 
previous studies. Nevertheless, as a result of the analysis 
of this amplitude equation, we also obtain a competi- 
tion between harmonic and subharmonic structures that 
agree in the limit of small modulation amplitudes very 
well with the results of the full numerical analysis of the 
Lengyel-Epstein model. 

This work is organized as follows. In Sec.|n]we present 
the Lengyel-Epstein model for a spatially modulated illu- 
mination and in Sec. lIIII we determine its stationary basic 
states and study their stability against small perturba- 
tions for a uniform and spatially modulated illumination. 
Some results of numerical simulations of the full nonlin- 
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ear model equations are presented in Sec. IIVI Close to 
threshold, in the so-called weakly nonlinear regime, the 
dynamical behavior of the Lengyel-Epstein model can be 
described in terms of an amplitude equation as discussed 
in Sec.E] The effect of the modulation on the threshold 
of this amplitude equation is investigated in Sec. IV Al bv 
using two different approaches given by a perturbation 
calculation and by a fully numerical solution of the gen- 
eral linear problem. The results are discussed in detail 
in Sec. IV HI where we also make a comparison with the 
thresholds obtained from a direct solution of the Lengyel- 
Epstein model. The work is finished with a summary and 
some concluding remarks in Sec. IVfl A detailed deriva- 
tion of the amplitude equation from the Lengyel-Epstein 
model is given in the Appendix. 



II. THE LENGYEL-EPSTEIN MODEL 

The starting point of our investigations on the effects 
of a spatially periodic modulated control parameter on a 
chemical reaction is the Lengyel-Epstein model 0, . 
This model describes two different instabilities of a spa- 
tially homogeneous chemical reaction, being either a Tur- 
ing instability to a stationary and spatially periodic pat- 
tern or a Hopf bifurcation to a spatially homogeneous but 
temporally oscillating reaction. Here we will focus on the 
Hopf bifurcation, which is preferred for similar diffusiv- 
ities of the reacting substances, and on the effects of a 
spatial modulated illumination in one spatial dimension. 
For this purpose, the model for the two dimensionlcss 
concentrations u{x,t) and v(x,t) is extended by a term 
describing a spatially varying illumination (j)(x), similar 
as in Ref. 331: 



dtu = a — cu — 4 
d t v = a \ cu 



1 + u 2 

uv 



1 + u 2 



-(f> + dlu, (la) 
+ <j> + ddlv\. (lb) 



The constants a, c, a and d denote dimensionless param- 
eters of the reaction diffusion model and the effect of 
an external illumination is introduced through the field 

fa), 



(a) = (f>o + M(x) 



(2) 



which can be identified as the control parameter of the 
system that is composed of a spatially homogeneous con- 
tribution cj)Q and a spatially varying part M(x). M(x) 
breaks the translational symmetry of the system, and we 
assume for reasons of simplicity a spatially periodic mod- 
ulation as described by 



M(x) = 2Gcos(2kx), 



(3) 



with the modulation amplitude 2G and the modulation 
wave number 2k. 



With the two vectors 



the matrix 



(4) 



C := 



d t + c-d 2 x 

-ere d t — add 2 ' 



and the nonlinear vector 



N := 



uv ( — 4 



1 + u 2 \—o 



(5) 



(6) 



a compact formulation of the two basic equations QJ be- 
comes possible, 



£u = N(u) + V, 



(7) 



which is especially useful for the amplitude expansion as 
outlined in the Appendix. 



III. BASIC STATES AND THEIR STABILITY 

The stationary basic state of the Lengyel-Epstein 
model is determined for a uniform illumination in 
Sec. IIII Al and for a spatially periodic illumination in 
Sec. IIII Bl In both cases we also investigate its stability 
with respect to a bifurcation to an oscillating chemical 
reaction. 



A. The spatially homogeneous case M(x) — 

In the case of a homogeneous illumination, i.e. for 
G = 0, the stationary solution of the Eqs. is given by 



Uq 



a — 50o 
5c 



a (1 



5w 



(8) 



It becomes unstable with respect to small oscillating per- 
turbations for an illumination strength (f> below a crit- 
ical value 0oc which is determined by a linear stability 
analysis. 

For this analysis we start with a superposition of the 
stationary, homogeneous basic state Uo and an infinites- 
imal perturbation Ui(x,t), 



u 



U + Ui = 



ui(x,t) 
vi(x,t) 



(9) 



as a solution of Eq. J2J , which is then linearized with re- 
spect to u\ and v\. The resulting two coupled differential 
equations 



£ui = A^oUi 



(10) 
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have the constant coefficient matrix 
M 



-4Ci 
-aCi 



4C 2 
aC 2 



with the matrix elements 
v (1 - u 2 ) 



C-2 



u 



1 + 



(11) 



(12) 



Equation IjlOjl may be further reduced by a mode ansatz 
of the form 



U| = A[ .k ] ' 



\t-\-iqx 



C.C. , 



(13) 



where c.c. denotes the complex conjugate and Eq de- 
scribes the ratio between the amplitudes of the two per- 
turbations ui and v\. The resulting two coupled linear 
equations have only solutions for a non-vanishing ampli- 
tude A 7^ 0, if the solubility condition 



\ + c + q 2 + 4d 4C 2 
1 -ac + ad X + adq 2 -\ 



aC 2 



= (14) 



is fulfilled. This condition determines the two eigenvalues 
X(q) as functions of the wave number q 



A± = ±^yj[q 2 (ad + 1) + 4Ci + c + aC 2 f - 4% 2 ) 



~ [ g 2 (ad+l) + 4Ci+c + (7C 2 ] 



(15) 



with 



h(q 2 ) = adq 4 + q 2 [ad(c + 4Ci) + crC* 2 ] + 5crcC* 2 . 

For a positive growth rate Re(A) and a finite imaginary 
part Im(A) ^ 0, the basic state Uo becomes unstable 
with respect to a Hopf bifurcation. The neutral curve of 
the wave-number dependent illumination strength 4>q (q) , 
which separates the stable from the unstable parameter 
range, is determined by the neutral stability condition 
Kc[X(q)] = 0. It is shown together with the correspond- 
ing Hopf frequency wo(?) f° r different values of the pa- 
rameters c in Fig. n Since a strong illumination of the 
chemical reaction suppresses the instability, the homo- 
geneous basic state uo is unstable for a given value of c 
within the area enclosed by the respective line in part (a) 
of Fig. ^ The maximum of each neutral curve 4>o(q) is 
given at q = that determines the critical illumination 
strength </>o c , below which the chemical reaction becomes 
oscillatory. In this case the eigenvalues in Eq. (|15|) may 
be further simplified to 

A± = ~ [4Ci +c + crC 2 ] 



±^y[4,C 1 +c + aC 2 ] -20acC 2 . (16) 

Since the parameters c, C 2 and a are all positive, also 
the product o~cC 2 is always positive and, therefore, the 
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FIG. 1: The neutral curves 4>o(q) in part (a) corre- 
spond starting from the top to increasing values of c = 
0.55, 0.65, 0.75, 0, 85, 0.95. In part (b) the curves for the Hopf 
frequency u>o(q) along <j>o(q) correspond starting from the bot- 
tom to increasing values of c. The remaining parameters are 
a= 12,<7 = 5,d = 0.8. 



eigenvalues given by Eq. i|16|) ar e either real with the 
same sign or complex conjugate. The latter case occurs 
if the condition 



c + 



4«o(l 



au 



(l + «§) 



2\ 2 



1 



20acuo 

1 +«n 



< 



(17) 



is fulfilled and the stability of the ground state Uo is then 
determined by Re(A±) = T(a,c,a,4>o) = —{AC\ +c + 
aC 2 )/2. The neutral stability condition r(a, c, a, 0o) = 
for the Hopf bifurcation is then in its explicit form given 
by 







125^ 



25(a - ha)4>l c 
a(3a 2 - haa - 125c 2 ) , 



25(5c 2 



2aa)(j) 0c 
(18) 



from which the critical illumination (f>o c (a,c,a) may be 
determined. The Hopf frequency at this critical value is 
described by the expression 



= ±^/hacC 2 = ± 



\ 



cr(a - 50o c ) 



1 



5c 



(19) 



Both, the critical illumination c/)q c and the Hopf fre- 
quency uj c , are plotted in Fig.[21as functions of the param- 
eter c and for three different values of a. With increasing 
values of c, the critical illumination decreases continu- 
ously up to the point 4>q c — 0. For c > c((f>Q C — 0) the 
stationary and homogeneous chemical reaction described 
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FIG. 2: In part (a) the threshold of the Hopf bifurcation, 
4>0c — 4>o(q — 0), and in part (b) the critical frequency u c = 
Ldo{q — 0) are shown as a function of the parameter c and 
for three different values of a. At small values of c, the Hopf 
bifurcation disappears and the homogeneous state uo becomes 
unstable with respect to a Turing instability. The data are 
determined for a = 5, d = 1. 



by Uo is always stable with respect to small perturba- 
tions. At small values of c, the Hopf bifurcation disap- 
pears and the basic state becomes unstable with respect 
to a Turing instability. 



B. Basic state in the presence of M(x) 

In order to determine in case of a spatially periodic 
illumination <f>{x) = <f>o + M(x) the stationary basic state 
do = (uq,vo) of Eq. 0, we use its time-independent part 
in the following form 

= a - 5cu - 5<p(x) + d 2 u - 4dd 2 v, (20a) 
= <z(l + u 2 ) - 5uv + (1 + v 2 ) (d 2 x u + dd 2 x v) . (20b) 

These inhomogeneously forced differential equations may 
be solved by the following Fourier ansatz for the two fields 
uq(x) and vq(x): 



m 



M 



u 



E 



il'2!:.: 



E 



H2!:.: 



(21) 



l=-M 



l=-M 



with an appropriate number M and the Fourier am- 
plitudes Ui and Vi of the expansion. The magni- 
tude of these amplitudes for I ^ is essentially de- 
termined by the forcing amplitude G. Since uq and 
Vq are real functions, we assume real amplitudes with 
Ui = U-i and Vi = V-i, respectively. Substituting the 
ansatz l|21f) into Eqs. I|20|) and, after projecting the equa- 
tions onto J dx exp (—ij2kx) in order to eliminate the 
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FIG. 3: The basic solutions uo (solid line) and vo (dotted 
line) of the Lengyel-Epstein model for a spatially modulated 
illumination rate cj>(x) — 4>o + 2G cos (2kx) with wave number 
k = 0.2 are shown in part (a) for G — 0.04 and in part (b) for 
G = 0.4. Parameters are a = 12, c = 0.55, d = 0.8, <j>o = 1-5. 



x-dependence, we obtain a set of coupled algebraic equa- 
tions for the determination of the unknown Fourier am- 
plitudes: 







= 



(a — 5^>o ) <5j,o — 5cllj — 5G (6j t i + <S 3 ',-i) 
-{2jk) 2 U ] + Ad{2jk) 2 Vj , (22a) 
a5j, +aJ2UiU j . l - 5^ U % Vj-i - (2jk) 2 U, 



^d(2 J k) 2 v J -j2J2 UiUmU - 



j — l — m 



[2mkf 



I m 



(22b) 



l m 



with j = —M . . . M. All sums in Eq. (|22bl) run from 
—M to M and the system of nonlinear equations in 
the amplitudes Uj and Vj can be solved by standard 
numerical methods. The basic spatially dependent so- 
lutions are then evaluated via Eq. (|2ip. One should 
note that according to Eq. I|22a|l . the Fourier amplitude 
Uq = (a — 50o)/(5c) corresponding to the spatially ho- 
mogeneous contribution to uq(x) is not changed by the 
forcing, cf. Eq. ©. 

For a given set of parameters, the two solutions uo and 
vq as given by Eq. (|21() are plotted in Fig. Ola) for the 
modulation amplitude G = 0.04 and in part (b) for G = 
0.4 while the modulation wave number is given by k — 
0.2. The field uq(x) is pictured as a solid line and Vo(x) as 
a dotted line. The spatial profile of the solutions shown in 
part (a) is dominated by the wave number of the forcing 
M(x) ex cos(2fcx). For increasing forcing amplitudes G, 
the weights of the higher harmonic amplitudes in the 
expansions given in Eq. I|21|l are amplified and the basic 
state Uo becomes fairly anharmonic as illustrated in part 
(b). Note the different scales in part (a) and (b). 



C. Threshold of the Hopf bifurcation in the 
presence of M(x) 

The spatially periodic basic state of the Lengyel- 
Epstein model as described by Uq = (uq(x), vq(x)) be- 
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comes unstable against infinitesimal perturbations w = 
(u>i, w 2 ) below a critical illumination rate (f>Q < <fi 0c (G, k). 
In order to determine this critical value, the basic state 
is separated from the small perturbation by the ansatz 

u = u + w . (23) 

After linearizing the basic equation Q) with respect to 
w, one obtains the following equation of motion 



Cw = M Q w, 
with the coefficient matrix 



M = 

and the abbreviations 



-4Cj 
-0-C1 



-4C 2 
-oC 2 



Cx{x) = 



(1 



-,2^2 



C 2 (x) 



"0 



(24) 



(25) 



(26) 



Equation (|24|) has formally the same form as for the ho- 
mogeneous case given by Eq. ljTU|) . Since Eq. (|2*1|) has 

spatially varying coefficients with a periodicity given by 
the forcing wave number, 2fc, the following Floquet-type 
ansatz for the small perturbations may be chosen with a 
complex parameter A: 



w 



N 

E 

l=-N 



e ilkx + c.c. 



(27) 



Substituting the ansatz l|27|l into Eq. I|24ll and using ad- 
ditionally the Fourier representation of C\ (x) and G 2 (x) , 



C x = 



A I 

E 

l = -AI 



^-.(1) 2ilkx 

O; e 



A I 



= E c > 



(2) lilkx 

I e 



(28) 



l=-M 



all terms can be sorted with respect to the linearly in- 
dependent exponential functions. To transfer the linear 
equation (|24|) into an eigenvalue problem for the constant 
coefficients Fi and Hi, one has to eliminate the remain- 
ing dependence on x by projecting the equations onto 
/ dxe~ % i kx . One finally ends up with the following sys- 
tem of equations: 



XFi = 



M 



A I 



j-2l 



4 E C\ 2) H^ 
(29a) 



l=-M 



<jk) 2 F 3 , 



l=-M 



A I 



£ C^H^-dijkfH.+cF, 



l = -M 
A I 



E c, 

l=-M 



(1) 



F, 



j-2l 



(29b) 



with j = —N . . . N. These equations can be written in a 
more compact form as two coupled sets of equations 

AF = AF+XiH, (30a) 



AH = A 2 U+I 2 F. 



(30b) 



where Ai and X, (i — 1,2) denote matrices of the di- 
mension (2N -I- 1) x (2iV + 1) and F and H include the 
(2N + 1) Fourier amplitudes of w\ and W2, respectively. 
Additionally, Eqs. (|30|l can be formally rewritten as an 
eigenvalue problem 



with 



C* = A* 



and C := 



Ax lx 

I2 A- 



(31) 



(32) 



The matrix C has the dimension (47V + 2) x (47V + 2). 
From Eqs. I|29|l one recognizes that the even and odd in- 
dices j are actually decoupled giving rise to two indepen- 
dent thresholds. These are the harmonic threshold 4>q c 
corresponding to harmonic perturbations w^x + ir/k) = 
w i( x ) (* — 1, 2) with respect to the forcing M(x) and 
the subharmonic threshold 0gJ corresponding to subhar- 
monic perturbations wf h (x + ir/k) — —wf (x). The 
larger one of these two thresholds determines whether 
spatially harmonic or subharmonic patterns emerge from 
the basic state iio via a Hopf bifurcation. 

Technically we solve the threshold problem as follows: 

(i) We keep besides the illumination 0o a h the other 
parameters a, c, d, a, G, k fixed and determine the basic 
states u and -Oo by solving Eqs. (|2*2*|l . 

(ii) Since the functions C\(x) and C 2 (x) depend accord- 
ing to Eq. i'2(j\) nonlinear ly on the basic state Uo, they 
are numerically evaluated on a discrete lattice and the 
Fourier amplitudes and C[ 2 \ as required in Eq. ij3"T|) . 
are obtained by a numerical Fourier transformation. 

(iii) The eigenvalue spectrum of the matrix C determines 
whether the basic state Uo is stable, i.e. if Re(A) < for 
all eigenvalues, or unstable, i.e. if Re(A) > at least for 
one eigenvalue A. Here A denotes the eigenvalue with the 
largest growth rate. The illumination (j>o is varied until 
the neutral stability condition Re(A) = is satisfied and 
this specific value of <fio determines the critical illumina- 
tion 4>q c , while Im(A) = ±lu c gives the Hopf frequency at 
threshold. 

Some results for the harmonic and subharmonic 
threshold as well as for the corresponding Hopf frequency 
are presented in Fig. 0] as functions of the amplitude G 
and in Fig. [3] as functions of the wave number k. Here, 
4>q c , are pictured as solid lines and <j)f^ , uj^ h as dashed 
lines. The harmonic solutions have according to Fig.QJa) 
the higher threshold for small forcing amplitudes G and 
from the basic state, a Hopf bifurcation occurs which is 
spatially harmonic with respect to the external modula- 
tion. However, the harmonic threshold drops below the 
subharmonic one for all values G > 0.077 and the bifur- 
cation from the basic state is changed to a spatially sub- 
harmonic pattern. The upper envelope of both threshold 
curves is the instability border below which the basic 
state Uo becomes unstable against small oscillating per- 
turbations. Subharmonic patterns are expected to occur 
in numerical simulations with random initial solutions 
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FIG. 4: In part (a) the threshold <^>o c of the harmonic insta- 
bility (solid line) and the threshold 0§c °f the subharmonic 
one (dashed line) are shown as functions of the forcing am- 
plitude G. Spatially subharmonic patterns are preferred in a 
wide range of G with cf>Q^ > (f>ij c . The dotted line indicates 
the lower end of the existence region of subharmonic patterns. 
The inset (b) displays the behavior of the thresholds for small 
values of G in order to enlarge the crossing region of the two 
curves 0o c (G) and 0qJ? (G). The two crosses in the inset mark 
the points at which some results of numerical simulations of 
Eqs. Q will be presented (see Fig. |5J. The Hopf frequency 
along the two threshold curves is shown in part (c) . Parame- 
ters are a = 12, c = 0.55, a = 5, d = 0.8, k = 0.2. 



within the area enclosed by the subharmonic threshold 
and by the dotted line in part (a) of the figure. Part (b) 
displays the two thresholds in the range of small values 
of G in order to magnify the neighborhood of the inter- 
section between (j)Q C (G) and 0g£(G). The two crosses in 
part (b) indicate the parameter values of (f>o and G, for 
which some nonlinear solutions of the Lengyel-Epstein 
model Q will be presented in Fig. [fj] The Hopf fre- 
quency along the two thresholds in part (a) is depicted 
in part (c). As can be seen from Fig.|S{a), the thresholds 
for the harmonic and subharmonic solution differ only 
slightly for small forcing wave numbers k and in the lim- 
iting case k — > 0, where the critical illumination is given 
by 4>oc = 4>0c(G — 0) + 2G, both thresholds coincide. 
On the other hand, for large forcing wave numbers where 
the modulation wavelength becomes smaller than the dif- 
fusion length, the system averages over the fast spatial 
oscillations ~ 1/fc, and the threshold approaches its un- 
modulated value 4>oc(G = 0) from below. The harmonic 
threshold <j>o c (k) has a pronounced minimum at k ~ 0.4. 
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FIG. 5: The thresholds for the harmonic instability 4>q c (solid 
line) and the subharmonic one (dashed line) are shown 
as functions of the forcing wave number k for a modulation 
strength G = 0.3. In a small range of k the subharmonic 
solution is preferred beyond the instability, <pQ C > <p 0c . In the 
limit k — > oo the harmonic threshold converges to 4>oc(G = 
0) = 1.47. In the opposite limit, i.e. k — * 0, both thresholds 
coincide and the critical illumination is given by (f>oc(G = 
0) + 2G. The Hopf frequency is shown in part (b). The 
remaining parameters are the same as used in Fig. [I] 

Here the bifurcation is almost suppressed, i.e. only a very 
weak illumination forces an instability of the stationary 
state Uo- For long- wavelength modulations k <C 1, the 
critical illumination is given by (f>o c > 4>0c{G — 0) and, 
therefore, the Hopf bifurcation already occurs in a range 
of the illumination </>o, where it does not appear with- 
out modulation. Again, in the /c-range where the har- 
monic threshold drops below the subharmonic one, the 
stationary basic state becomes unstable against spatially 
subharmonic perturbations. The corresponding Hopf fre- 
quency is displayed in part (b). 



IV. NONLINEAR SOLUTIONS 

The determination of the time evolution of the non- 
linear solutions of Eqs. Q for spatially periodic bound- 
ary conditions is performed by a pseudospectral method. 
From a numerical point of view, it has been proven useful 
to separate the stationary basic state iio from the oscilla- 
tory contribution w in order to simulate the equations of 
motion. After inserting the ansatz \TM in the Lengyel- 
Epstein model Q, one obtains the following governing 
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FIG. 6: The left column shows harmonic solutions for the 
fields w\ (top) and v = vq + w\ (bottom) and the right 
column subharmonic solutions for the same fields, i.e. wf 
(top) and v = i>o + wf (bottom). Solid lines indicate 
the spatial modulation M(x). The illumination is given by 
4>o — 1.42 and the forcing amplitude by G — 0.025 in part 
(a) and by G = 0.15 in part (b). Further parameters are 
a — 12, c = 0.55, a = 5, d = 0.8, k = 0.2. The parameter 
values in the <j>o — G plane are marked by crosses in Fig. Efb). 

equations for the two fields w\ and w 2 : 

d t wi = —cwi + d 2 x w\ 

\(uq + w 1 )(v + w 2 ) u v 1 v 

-4 — 77 75 79 , (33a) 

[ l + (u +wi) 2 l + «oJ 

dtW2 = <j(cwi + dd^w 2 ) 

{u + wi)(v + w 2 ) uqV 1 foo ,, 

—a ; — K .(oob) 

l + (lt + Wl) 2 l + "oJ 

The effect of the modulation M(x) enters these equations 
via the basic states uq and vo- A great deal of simulations 
with random initial conditions were performed in order 
to verify the onset of the Hopf bifurcation numerically. 
The results are in excellent agreement with the threshold 
curves for 0q c (G) and 0gJ (G), respectively, in Fig.^Ja). 

Two types of nonlinear solutions obtained by numer- 
ical simulations are presented in Fig. |SJ where the left 
column shows the time evolution of a spatially harmonic 
solution for the field w 2 (top) as well as for the super- 
position v — vq + w 2 (bottom) occurring at <fio — 1.42 
and G = 0.025. The same fields are shown in the 
right column, i.e. w 2 h (top) and v = vci + w 2 h (bot- 
tom), in the case of subharmonic solutions occurring at 
G = 0.15. The parameters for these simulations in the 
(f>o — G plane are marked by the crosses in Fig. 01b). 
The solution w 2 (x,t) describes spatially modulated os- 
cillations whose spatial profile becomes more pronounced 
when the basic solution vq is included. The time evolu- 
tion of the field w 2 h resembles that of a standing wave 
with twice the wavelength of the forcing and, indeed, it is 



well described by a superposition of three standing waves 
wf(x, t) = Y% =1 Bj sin (ujft + if,) sin [(2j - l)kx] with 
real amplitudes Bj and phases ifj. On the bottom right 
the full solution v = io + w 2 h of the Lengyel-Epstein 
model is shown. Note, the basic solutions vq, uq have the 
same periodicity as the spatial modulation. The nonlin- 
ear solutions for the other fields wi and u = iio + w\ look 
very much like the ones shown in Fig. and they are 
therefore not presented here. Starting the simulations 
with spatially subharmonic states and decreasing the il- 
lumination cf>o continuously, the subharmonic pattern is 
stable against small perturbations up to the dotted line 
in Fig. 2Ja). On exceeding this border line, the subhar- 
monic solution becomes unstable and spatially harmonic 
patterns emerge. 

Close to the threshold of the Hopf bifurcation, where 
the amplitude of the oscillations is small, the dynam- 
ics of the Lengyel-Epstein model may be described by 
a so-called amplitude equation as discussed in the next 
section. 



V. AMPLITUDE EQUATION 

Below the critical value (f>o c of the control parameter 
0o, the basic state Uo = (uq,vo) of the Lengyel-Epstein 
model Q becomes linearly unstable against small oscil- 
latory perturbations Ui as described in Sec. IIII Al The 
magnitude of Ui beyond threshold is restricted by nonlin- 
ear terms in Ui which are of cubic order for a supercritical 
bifurcation. In this case one may derive a universal am- 
plitude equation for the slowly varying amplitude A(x,t) 
in order to describe the dynamics close to the threshold 
[38| . The technique for this derivation is a multiple scale 
analysis where the full solution Ui is decomposed into a 
fast varying oscillation oc exp (iuj c t) with the frequency 
w c and a spatially and temporally slowly varying ampli- 
tude A(x,t): Ui = A(x, t)u e exp (iu} c t) + c.c. 

The usage of amplitude equations is a well established 
method to characterize the universal properties of a pat- 
tern at small amplitudes close to its threshold. A par- 
ticularly well-known amplitude equation is found for a 
spatially homogeneous Hopf bifurcation, which has been 
investigated very intensively over the recent decades and 
a recent review of this subject is given by Ref. [3^ . 

For the derivation of the amplitude equation one in- 
troduces as an expansion parameter the small distance 
to the threshold e = and the expansion holds in 

the range A cx e 1 / 2 . Here we assume additionally that 
the modulation M(x) is also of the order e. For this case 
we generalize the amplitude equation for an oscillatory 
bifurcation by including the spatial modulation M(x): 

r d t A = e{l + i Cl )A + il{l + ib)dlA (34) 
+si (1 + is 2 ) MA- g (1 + ic 2 ) | A | 2 A . 

An explicit derivation of this equation from the basic 
Eqs. is given in the Appendix and a special case of it 
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is given in Ref. All coefficients in Eq. I|34[l describe 
physical quantities, such as the relaxation time To, the 
linear and nonlinear frequency shift c\ and C2, respec- 
tively, the coherence length £0 and the linear frequency 
dispersion £q6. For g > the bifurcation is supercritical 
and otherwise subcritical. The analytical expressions for 
all these coefficients in terms of the parameters of the ba- 
sic equations JTJ are rather lengthy and, instead of giving 
their analytical forms, we have plotted them in Fig. as 
functions of the parameter c and for three different values 
of a. 

It is worthwhile mentioning that apart from the coef- 
ficients si and S2, all the other linear coefficients can be 
calculated from the dispersion relation \((f)o , q 2 , . . .) — 
Re(A) ± ilm(A) of the Lengyel-Epstein model given in 
Eq. JTSJ by the following expressions [H, Ejj 



1 



TO 



e 



(j> 0c dRe(X)/dcl>o 
1 d 2 6n 



ci 



duj 
ocpo 



20oc dq 2 



b=- 



T d 2 LO 



(35a) 
(35b) 



Here u> = Im(A) and the derivatives are evaluated at 
the critical values 4>o c ,q c ,uj c . For a vanishing modula- 
tion these coefficients describe the linear properties of 
the Lengyel-Epstein model near threshold. It is however 
indispensable to carry out the perturbation expansion in 
order to determine the linear coefficients S\ and S2 as well 
as the nonlinear coefficients g and C2 as functions of the 
parameters of the basic equations. 

The term iec\A in Eq. (|34|l can be removed by the 
transformation A = e~ l£Clt A. For convenience we can 
scale out the coefficients tq,^q, si, g by a suitable choice 
of time, space and amplitude scales 



t' = t/r , 
^'= 5 1/2 A, 



x 
G' 



x/£,o , 
siG, 



(36) 



and one obtains the following rescaled amplitude equa- 
tion: 



d t A = eA + (1 + ib) d%A + (1 + is 2 ) MA 
-(1 + ica) I A\ 2 A, 



(37) 



where we have kept for simplicity the same symbols 
for the scaled quantities. The amplitude equation is 
invariant under an arbitrary phase transformation as 
A — > Aexp (itp). 



A. Determination of the threshold 

We investigate in this section how the spatial modula- 
tion M(x) changes the bifurcation scenario from the basic 
state A = of Eq. (|37|l into a spatial pattern. In the ab- 
sence of the modulation, i.e. G — 0, the linear part of 
Eq. is solved by A = Fe xt+iqx and this ansatz leads 
in the neutrally stable case Re(A) = to an expression 




FIG. 7: The coefficients of the amplitude equation il' J > H are 
plotted as functions of the parameter c of the Lengyel-Epstein 
model and for three different values a: a = 14 (solid lines), 
a — 12 (dashed lines) and a = 10 (dotted lines). The critical 
illumination <j>0c and the Hopf frequency u> c along these curves 
are shown in Fig. [5] Further parameters are a = 5, d = 1. 



for the neutral curve £o(q) and the frequency dispersion 
uj (q) = Im(A), 



vo(q) = ~bq 2 



(38) 



Minimizing eo(q) determines the threshold e c = 0, the 
critical wave number q c = and the frequency oj c = 0. 
In the presence of M(x) the bifurcation properties of 
Eq. l|3T|) are changed, as illustrated in the following by 
a perturbation calculation for small modulation ampli- 
tudes. 



1. Perturbation method for small amplitudes of M(x) 

For small forcing amplitudes G « 1 we introduce a 
small expansion parameter 77 1 with M(x) = rjM(x). 
Since the amplitude equation (|37|l is of first order with 
respect to time, the solution of its linear part has an 
exponential time dependence and for small values of the 
modulation amplitude, G — rjG, the linear solution may 
be expanded in powers of the modulation strength r\ 



A = e xt (A + 77A1 + j] 2 A 2 + ■■■) 



(39) 



Applying the neutral stability condition Rc(A) = 0, the 
perturbation in Eq. 139|) does neither grow nor decay, 
which separates the parameter regime where the basic 
state A = is stable from the range where A = is 
unstable. 

The threshold is shifted due to the modulation and 
therefore, the control parameter e c and the frequency 
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uj c = Ini(A) are also expanded with respect to r] 



and for the subharmonic case by 



e c = 4°) +?7 eW+r, 2 4 2) 



(40a) 
(40b) 



The expansions given in Eqs. I|39l) and l|40|) provide the 
following hierarchy of equations defining the neutral sta- 
bility of the basic state A = 0: 



77° : C A = , 



(41a) 



77 1 : £ Ai = [e^ +1^ + {l+is 2 )M]A , (41b) 



rf : C A 2 = [ei 1) +w£> + {l + is 2 )M]A 



(i) 



+ [ S V> +iu,W]A 



(2)1 



(41c) 



with the linear operator £ = d t — iu>c — £c°' ) — (l + ib)d 2 . 
These equations may be solved by a spatial dependence 
which is either harmonic 



A a = F , 

A x = F 2 e 2ikx + F_ 2 e- 2ikx 



(42a) 
(42b) 



or subharmonic 



A 1 



F_ie" 



Fie 

? 3 



(42c) 
(42d) 



with respect to M(x). In order to distinguish between the 
harmonic and subharmonic results we introduce Shi^h 
for the harmonic case and E s h^sh for the subharmonic 
one. From a solubility condition on the right hand side 
of Eqs. (|41b|> and (|41c|> the corrections to e c and cu c may 
be calculated. The solution A\ of Eq. i|41b[) is given for 
the harmonic ansatz by 



G(l + is 2 )F 



4k 2 (l + ib) [ + h [ ' 

and in the case of the subharmonic ansatz by 



Ar 



G(l + is 2 ) , 3tkx 



8k 2 (l+ib) 



F ie MKX + F_ie 



— 3ikx\ 



(44) 



whereas the solution A 2 of Eq. (|41c|) is not needed ex- 
plicitly to determine the corrections and u4 , re- 
spectively. The expansions of the threshold e c and the 
frequency ui c up to order 0(rf) are given for the harmonic 
case by 



,G 2 (l-s 2 + 2bs 2 ) 
2fc 2 (l + 6 2 ) 
2 G 2 [b{l-s 2 )~2s 2 ] 
2fc 2 (l + 6 2 ) 



(45a) 
(45b) 



, 2 n 2 G 2 {l-s 2 + 2bs 2 
e sh = k -r)G-r) 



u s h = bk 2 - f]s 2 G + T) 



8fc 2 (l + 6 2 ) 
2 G 2 [6(1 - s 2 ) - 2s 2 ] 



8fc 2 (l + fe 2 ) 



(46a) 
(46b) 



If one increases the control parameter in Eq. I|37[l from 
below, the lowest of the two thresholds Eh and e s h deter- 
mines whether the basic state A = is unstable against 
spatially harmonic solutions, i.e. if £h < £ s h, or spatially 
subharmonic solutions, i.e. if e s h < £n- By replacing 
k — > £ok and G — > siG in the expressions l|45|l and l|46|l . 
the thresholds £h,sh and the frequencies tOh,sh for the am- 
plitude equation (|34|l follow. 



2. Numerical determination of the threshold 

Because of the periodically varying coefficient M(x) in 
Eq. I|37[) . the general linear solution may be represented, 
similar as in Sec. IIII CI by a Floquet-Bloch expansion 



A(x,t) 



N 

\t+iqx \ * 
j=-N 



F 3 e i32kx (with < q < 2k) . 



up to an appropriate number N . Without spatial forc- 
ing all coefficients but Fq vanish. The perturbation in 
Eq. (|47|l grows for a chosen parameter combination if 
Re (A) > and it decays if Re(A) < 0. We are again pri- 
marily interested in the neutrally stable case Re(A) = 0, 
separating the stable from the unstable regime. After in- 
serting the ansatz for A into the linear part of Eq. 13711 , 
the explicit ^-dependence is removed by multiplying the 
equation with exp (—U2kx) (I = —N, . . . , N) and inte- 
grating with respect to x. From this procedure the eigen- 
value problem 



AF = pF 



F = {F_ N , 



,F N ) 



(48) 



follows, where the matrix A is a band matrix of width 
(27V + 1) with the coefficients 

Ai,i = {I + tb){q + 2lk) 2 , (49a) 
Ai,i- 2 = -(l + is 2 )G, (49b) 
A M+ 2 = A U - 2 . (49c) 

From a solvability condition for the homogeneous system 
of equations (|4*5|l , 



det(„4 - pi) = f(G, k,q...)=0. 



(50) 



(X is the unity matrix) the eigenvalues pi are deter- 
mined as functions of the parameters and they can be 
sorted in ascending order with respect to their real parts 
Re(pi). Keeping all parameters besides q fixed, the neu- 
tral curve £o(q) and the frequency uJo(q) can be de- 
termined from the eigenvalue with the lowest real part 
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p = min[Re(pi)]. Minimizing eo(q) gives then the thresh- 
old e c , the critical wave number q c and the critical fre- 
quency uj c = Wo(q — q c )- In this way, we found that the 
minimum of the neutral curve is either given at q c = 
determining the harmonic threshold £h, or at q c = k 
determining the subharmonic threshold e s h- The har- 
monic and subharmonic contributions of the ansatz in 
Eq. H47fl with respect to M(x) separate for the linear 
part of Eq. (|37() , and the two thresholds Eh and e s h may 
be calculated independently. 

In order to obtain the analogous eigenvalue equation 
with respect to the unsealed amplitude equation given 
in Eq. I|34l) . one has to replace in Eqs. (|49|l the wave 
numbers q and k by £09 and £,ok, respectively, as well as 
the amplitude G by s±G. 



B. Results 

One interesting question is the location of the boarder 
separating the parameter range where the harmonic pat- 
tern is preferred at threshold from that range where the 
subharmonic pattern is favored, whereby the borderline 
is determined by the condition Eh — £ s h ■ In terms of our 
perturbational results as given in Sec. IV A II this condi- 
tion leads to a second order polynomial in the modulation 
amplitude G with its two solutions 



G± = 



4fc 2 (l + fr 2 ) 
3(1 - si + 2bs 2 ] 



± 



2k 2 (l 



3(1 - s 2 2 + 2bs 2 ) 



. 1 - si + 2bs 2 



l + b 2 



(51) 



There is a finite range in G where subharmonic solutions 
are preferred, namely if the following inequality 



4 (1 + b 2 ) - 6 (1 - 4 + 2bs 2 ) > 



(52) 



is fulfilled. The harmonic threshold is the lowest one for 
modulation amplitudes G < G_ and G > G+. In the 
finite range G_ < G < G+ the subharmonic threshold 
e s h drops below the harmonic one. The two amplitudes 
G + and G_ according to the formula (|ST|) are plotted in 
Fig. [8] as functions of the coefficient s 2 ■ The ranges in 
which subharmonic or harmonic solutions are preferred 
are marked by s and h, respectively. 

In the parameter range s 2 > b + yl + b 2 and s 2 < 
b — vl + b 2 the harmonic threshold has a positive cur- 
vature as a function of G, d 2 £h/dG 2 > 0, which can be 
readily verified from Eq. (|45a|l . Consequently, at small 
values of G the threshold Eh is shifted upwards while 
the subharmonic threshold is dominated by the linear 
decrease e s h oc k 2 — G. These trends of the two thresh- 
olds apparently promote the appearance of subharmonic 
patterns. 

For arbitrary values of the modulation amplitude G 
and wave number k the linear stability of the basic state 
A = of Eq. I|37fl has to be determined numerically 



O 

C 

I 

O 



I 0.2 




coefficient s 2 

FIG. 8: The two solutions G+ (solid lines) and G_ (dot- 
ted lines), as given by Eq. I5H . are shown as a function of 
the coefficient S2 and for the parameters k — 0.25,6 = 1.5. 
The areas with pure harmonic and subharmonic solutions at 
threshold are marked by h and s, respectively. 



by solving the eigenvalue problem |(48j. The harmonic 
threshold Sh and the subharmonic threshold e s h are cal- 
culated separately and some results for them as well as 
for the corresponding critical frequencies cuh and Lu s h , re- 
spectively, are shown in Fig.[5]for two sets of parameters. 
The harmonic branches Sh , ^>h (solid lines) and subhar- 
monic branches e s h,w s h (dashed lines) have been calcu- 
lated for N — 32 Fourier modes, which has been proven 
to be a reasonable approximation. The symbols in part 

(a) indicate the results of the perturbation calculation as 
given in Eqs. (|45|l and l|46(l . which are in good agreement 
with the numerical results for small forcing amplitudes 
G. The points of intersection between the branches £h 
and s s h are roughly given by G_ ~ 0.068, G + ~ 1.53 
in part (a) and by G_ ~ 0.028, G+ ~ 0.12 in part 

(b) . For comparison, the perturbation calculation yields 
G- ~ 0.068, G+ ~ 0.68 for the parameters in part (a) 
and G_ ~ -0.048, G+ ~ 0.027 for the ones in part (b). 
As mentioned above, the sign of the curvature d 2 Eh/dG 2 
at small amplitudes G may be changed by varying the co- 
efficients b and s 2 , which can be recognized by comparing 
the course of Eh in part (a) given for (6, s 2 ) = (5, 0.5) and 
in part (b) given for (6, s 2 ) = (0.01, 3). 

Harmonic solutions are preferred for large modulation 
wave numbers fc, and in the limiting case k — > oo, the 
harmonic threshold approaches Sh — while the subhar- 
monic threshold diverges in this limit, being in agreement 
with the expressions given in Eqs. (|45a|l and (|46a|l . In 
the opposite limit, i.e. k — > 0, the thresholds Eh and s s h 
tend to Eh,sh = — 2G, which has to be computed numer- 
ically. We remind the reader that the thresholds <pQ C (k) 
and <Poc(k), as obtained for the Lengyel-Epstein model, 
exhibit qualitatively the same behavior as a function of 
k, cf. Fig. Efa). 

In Fig. the harmonic threshold Sh (solid line) and 
the subharmonic threshold e s h (dashed line) are plotted 
as functions of the linear coefficient s 2 (top) and b (bot- 
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amplitude G amplitude G 

FIG. 9: The harmonic (solid lines) and subharmonic (dashed lines) threshold as well as the corresponding frequencies are 
plotted as functions of the forcing amplitude G for two sets of parameters, in part (a) for k = 0.25, b — 5, S2 = 0.5 and in 
part (b) for k = 0.25, b — 0.01, S2 = 3. In order to generate the numerical results N — 32 modes in the Floquet expansion 
1471 1 have been used. Subharmonic patterns occur in the G-range where e s h < Sh- The second point of intersection between 
the harmonic and subharmonic branch is located at G ~ 1.53 in part (a) and at G ~ 0.12 in part (b). For comparison, the 
symbols in part (a) are obtained by the analytical expressions given in Eqs. I|45^ and 1460 . 



torn). For large values of the modulus \a%\ the subhar- 
monic threshold drops below the harmonic one and spa- 
tially subharmonic solutions appear at threshold of the 
Hopf bifurcation. Regions solely supporting harmonic 
or subharmonic patterns are found in the range s and 
h, respectively. The points of intersection between both 
thresholds are roughly given by S2 — ±2. For compar- 
ison, S2 — ±1.84 is obtained from the perturbation cal- 
culation. According to the bottom part of the figure, 
subharmonic solutions are preferred in a finite range of 
the parameter b similar as the G-range in Fig. [!|] Again, 
within the regions marked by s and h only subharmonic 
or harmonic patterns are found. 

From the inequality given in Eq. 1)52(1 one may easily 
deduce the following limiting cases. For 6 = subhar- 
monic solutions are favored in the range s\ > 1/3 and for 
S2 = in the range b 2 > 1/2, while b = S2 = seems to 
lead to a contradiction. Therefore, one of the two linear 
coefficients b and S2 is necessary for the occurrence of 
spatially subharmonic patterns. 



1. Comparison with the Lengy el-Epstein model 

In order to be able to compare directly the results for 
the harmonic and subharmonic thresholds with those ob- 
tained from the basic equations QJ, we have to con- 
sider the linear part of the unsealed amplitude equa- 
tion (|34|l . which reflects the typical length scale, £o> time 
scale, tq, and amplitude scale, Si, of the Lengyel-Epstein 



model. In this context one may raise the question under 
which conditions are the thresholds Eh,£ s h determined 
by the amplitude-equation approach, and the thresholds 
4>q c ,4>oc determined by the Lengyel-Epstein model, in 
(good) agreement? Such a comparison of the thresholds, 
e.g. as a function of the modulation amplitude G, gives 
also some insights into the validity range of the amplitude 
equation which is a priori unknown. 

For this comparison, we have to replace in the eigen- 
value problem in Eq. I|48|) the wave numbers q and k by 
£o9 and £ok as well as the amplitude G by siG. Further- 
more, the values of the linear coefficients tq, Ci, £oj b, Si 
and S2 are determined by the parameters of the Lengyel- 
Epstein model as shown, for instance, in Fig. [7| 

Figure^2 sn ows the harmonic threshold (triangles) and 
the subharmonic one (squares) as given for the Lengyel- 
Epstein model as well as the related thresholds obtained 
from the amplitude equation (lines). The latter ones are 
calculated with respect to the illumination rate by us- 
ing the formula 4>q c = 4>o c {G = 0)(1 — e c ). According 
to this figure, one finds a qualitatively similar behavior 
of the associated thresholds. For comparison, the rel- 
ative error between the harmonic thresholds is roughly 
given by 5% and between the subharmonic thresholds by 
6% calculated for G = 0.06. These deviations become 
smaller for decreasing values of the wave numbers k or 
by decreasing the modulation amplitude G. Moreover, 
the finite G-range in which the subharmonic solution has 
the higher threshold is also in good agreement with the 
range given by the amplitude approach. 
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FIG. 10: The dependence of the harmonic (solid line) and 
subharmonic (dashed line) threshold on the coefficient S2 
is shown in the top part for a set of parameters given by 
G = 0.04, k = 0.25, b = 0.01. The bottom part shows both 
thresholds as function of the coefficient b and for the param- 
eters G = 0.05, k = 0.24, S2 = 3. The range where harmonic 
solutions are preferred is marked by h and that of subhar- 
monic ones by s. 

Another interesting comparison between both models 
is the change of the sign of the curvature of the har- 
monic threshold, which has been found for the amplitude 
equation, cf. Eq. I|45a|) . This significant behavior of the 
threshold has also been observed for the Lengyel-Epstein 
model in a parameter regime determined by the ampli- 
tude equation. If the curvature of 4>o c {G) is positive the 
Hopf bifurcation also occurs for larger illuminations <j>o, 
being in contrast to the situation shown in Fig. EJa) or 
Fig. El 



VI. SUMMARY AND CONCLUSIONS 

We have investigated the effects of a spatially periodic 
modulated control parameter on an oscillating chemical 
reaction described by the Lengyel-Epstein model. This 
reaction exhibits a supercritical Hopf bifurcation and it 
provides due to its photosensitivity a simple approach 
to study the response of the reaction with respect to a 
spatially modulated illumination. 




0.66- 
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FIG. 11: Comparison between the harmonic (h) and sub- 
harmonic (s) thresholds, which have been calculated for the 
Lengyel-Epstein model (symbols) by solving the eigenvalue 
equation l-'i 1 > and for the amplitude equation (lines) by solv- 
ing the eigenvalue equation 1481 . For the latter one, we 
have replaced the wave numbers q and k by £o<Z and £ok 
as well as the amplitude G by siG. The modulation wave 
number is k — 0.05. The other parameters are given by 
a = 10, c = 0.55, a = 5, d = 1 leading to the coefficients of the 
amplitude equation 1341 as given in Fig. |7| (dotted line). 

We find that in the range of intermediate values of 
the modulation amplitude G the bifurcation to the os- 
cillatory chemical reaction is subharmonic with respect 
to the external modulation whereas for small and large 
modulation amplitudes the bifurcation is harmonic. Be- 
yond the bifurcation point the subharmonic solution is 
preferred in a finite range of the control parameter before 
a transition to the harmonic pattern takes place at larger 
values. The related results of the stability calculations of 
the basic state with respect to small perturbations are 
summarized in Figs. 01 and [S] 

Close to the threshold of the Hopf bifurcation, an am- 
plitude equation is presented which is an extension of the 
well-known complex Ginzburg-Landau equation for spa- 
tially homogeneous, oscillatory bifurcations. This equa- 
tion is generic for oscillatory systems near threshold un- 
derlying a spatially varying control parameter and it may 
therefore also be found for other systems having the same 
symmetry properties as the considered Lengyel-Epstein 
model. 

The stability limit of the basic state of the amplitude 
equation has been determined by a perturbation calcu- 
lation and by solving the general linear problem numer- 
ically. Good agreement of the two approaches is found 
for small forcing amplitudes. It has been shown that in- 
termediate forcing amplitudes lead also to subharmonic 
solutions, while weak and strong forcing amplitudes favor 
harmonic solutions. A rough estimation of the validity 
range of the amplitude equation has been given by com- 
paring the harmonic and subharmonic thresholds with 
those obtained for the Lengyel-Epstein model itself. We 
have found that the amplitude equation describes the 
linear properties of the Lengyel-Epstein model to a great 
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extent for long-wavelength modulations and small forcing 
amplitudes. 

According to our results we expect in experiments on 
chemical reactions, which are described by the Lengyel- 
Epstein model, also a transition to subharmonic patterns 
induced by a spatially periodic illumination. 

Instead of considering stationary forcing, the effect of 
spatiotemporal forcing on an oscillating chemical system 
may also be investigated. Here, the forcing has the form 
of a traveling wave similar as introduced recently to study 
the effects of spatiotemporal forcing on stationary Tur- 
ing patterns pa. l2Sj as well as on oblique stripe patterns 
in anisotropic systems [3(j. This special type of forcing 
breaks a further symmetry of the system, the reflection 
symmetry, which may induce a more complex spatiotem- 
poral behavior to which forthcoming work is devoted. 



APPENDIX: SCHEME FOR THE DERIVATION 
OF THE AMPLITUDE EQUATION 

For the derivation of the amplitude equation (|34|l . a 
small reduced control parameter is introduced, 



uqc - <Po 

4>0c 



(A.l) 



which is a measure for the distance from the threshold of 
the Hopf bifurcation. At threshold the linear solutions of 
the Lengyel-Epstein model may be written as 

ui = Au e e^ c * + c.c. , (A.2) 

where A describes the common amplitude of the two 
fields ui,vi and the eigenvector u e = (l,Eo) T describes 
their amplitude ratio, cf. Eq. (|13|) . 

For the perturbation expansion we introduce slow time 
and space variables [3£| 

X = e 1/2 x, T l= e 1/2 t, T = et (A.3) 

and the solution Ui may be written as a product of a 
slowly varying amplitude and a fast oscillating exponen- 
tial function 



ui = A(X, T, Ti)u e e lUct + c.c 



(A.4) 



According to the chain rule one may replace time and 
space derivatives by the following expressions 



a t dt + e 1/2 a_ 



d x ^d x + e 1 ' 2 d x . (A.5) 



Note that dt and d x on the right-hand side only act on the 
rapid dependences. We further assume small modulation 
amplitudes M(x) = eM(x) with M(x) cx 0(1) and wave 
numbers k = e 1 / 2 ^ with k cx 0(1). An attribute of 
the ± symmetry of a supercritical oscillatory bifurcation 
is the power law for the oscillation amplitude A ~ ^/i. 
Accordingly, the solutions of the basic equation Q are 
expanded near threshold with respect to powers of e 1 / 2 , 



u 



•eVa Ul 



£U 2 



e 3/2 u 3 



(A.6) 



The components of Uo describe the stationary solu- 
tions as given in Eq. iJSJ and the components of ii! = 
A\\i e exp (iu) c t) + c.c. provide the linear oscillatory solu- 
tions at threshold as discussed in Sec. IIII Al Note that 
the amplitude A\ and the amplitude A in Eq. I|A.4|) only 
differ by a factor of e 1 / 2 . The expansions for C, N and 
V in Eq. J7J are given by 

C = C + £ 1/2 A + eC 2 , (A.7a) 

N = N + L + £N 2 + £ 3/2 N 3 + --- , (A.7b) 

with L = e 1/2 MqVli + eM u 2 + e 3/2 M u 3 H , 

V = V + sV 2 , (A.7c) 

and with the explicit expressions of the linear operators 
£o,£i as well as L 2 . 



C 
Hi 
C 2 



0, 
d Tl 



-trc d t 



2d x d x 
dr, 

d T -d\ 

8t — crdd\ 





2add x dx 







(A. 



The nonlinear vectors in the expansions l|A.7b|) 
and ljA.7cfl are given by 



N 
N 2 

N 3 

7V 3 
and 



U V 



-4 



(Eiuivi + E 2 u\) , 

7V 3 , with 
[Ri{u\v 2 + u 2 vi) + R 2 uiu 2 + R3V111I + R^ul] 



V = 



a - <p 0( 



V 2 = 



where the abbreviations 
1 



E! = 

Ri = 
R* = 



(1+^) 3 



(1 



Ri 



90c 

-<J<f>0c 



V (ul 



M _ 



3w n 



(A.9) 



(1+^) 
-2u uo (3 



2u 2 



4) 



(1 



(A.10) 



(1 



Ra — 



6ulv a - U^Vg - V 
(l + «§) 4 



have been introduced. Inserting all expressions into 
Eq. |0 yields at successive orders of e 1 ! 2 

(A. 11a) 
(A.llb) 
(A.llc) 
(A.lld) 



e° : 


-C u 


= N + V , 


e 1 ' 2 : 


£ ui 


= M ui , 


e 1 : 


£ u 2 


= Mqu 2 - £im + N 2 + V 2 


£ 3 / 2 : 


C u 3 


= Movl-3 + N 3 - £ 2 ui , 
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Equation ijA.llafl determines the basic state as given by 
Eq. ||SJ| while Eq. (|A.llb|l reproduces the threshold condi- 
tion H14J with the solution known already from Eq. (|15|1 . 
Since (Co — Mq)ui = holds due to Eq. jA.llbjl . the 
contribution C\Ui on the right-hand side of Eq. (|A. 1 lc|) 
has to vanish leading to the constraint 8t ± Ai = 0. In 
order to solve the inhomogeneous equation i|A.llcfl . we 
choose the following ansatz: 



u 2 



A 
Bo 



A 2 \ 2iw c t 

B 2 



(A.12) 



After inserting Ui and u 2 into Eq. (jA.lldfl . there is no 
need to solve this equation explicitly. Projecting instead 
the whole equation onto lit, where u{ is the solution 
to the adjoint equation of l|A.llb[) . the inhomogeneity 
on the right-hand side of Eq. l)A.lld|) yields a solubility 
condition. By rescaling back to the original units of the 
space and time variables and also to the amplitudes M = 
eM and A — e^-^Ax, this solvability condition provides 
the final form of the amplitude equation for A: 



For the four unknown amplitudes we obtain 

0oc - M 



.4, 



B 9 



A 

c 

Cx(<p 0c -M) Ex (E + E$) + 2E 2 2 
Bo = — — | -All , 



C 2 c C 2 

-8iuj c (ExE + E 2 ) 

5<tcC 2 - 4a)2 + 2ico c (c + 4Ci + aC 2 

<r(5c + 2w e ) (ExE +E 2 ) 
5<jcC 2 - Aw 2 c + 2iiv c (c + 4Ci + aC 2 



A 2 



^.(A.13) 



To d t A = e{l+icx)A + £l{l + ib)dlA (A.14) 
+sx (1 + is 2 ) MA -g(l + ic 2 ) \A\ 2 A. 



All coefficients are given in terms of the parameters of the 
basic equations JQ) and we have plotted them in Fig. 
as functions of the parameter c and for three different 
values of o. 
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